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We consider the reaction yd— >ir°d in a wide energy range around and above the 77-meson pho- 
toproduction threshold at backward CM angles of the outgoing pion. Our theoretical analysis is 
essentially motivated by the recent measurements of the CLAS Collaboration at Jefferson Lab, 
where this kinematical region of the reaction has been thoroughly studied for the first time and a 
cusps in the energy dependence of the differential cross section in the region of E- t ~ 600 — 800 MeV 
has been observed. Our preliminary and qualitative analysis, based on single- and double-scattering 
diagrams, shows that the observed structure can be explained by the contribution of the double- 
scattering diagram with intermediate production of the rj meson. The effect, to a considerable 
extent, is formed due to the contribution of JV(1535) resonance to the amplitudes of subprocesses 
on the nucleons. 
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has recently been studied [lj at Jefferson Lab using the CLAS detector. The experiment was carried out in a wide 
range of photon energies E 1 — 0.5 — 2.0 GeV at large CM scattering angles of the outgoing pions. A new phenomenon 
was observed. 

The process of pion photoproduction on the nucleon jp—urN has been theoretically investigated over a long period 
of time (see, for example, Refs. Ill E3> El 13 ) • The reaction Q on the deuteron, has been also considered in a number 
of papers [f| Q. During the last years, the 77-photoproduction processes have also been actively studied (see, for 
example, Refs. 0,0 for yp^rjp, for yd^rjd, and 0] for 77 photoproduction on light nuclei). 

The study of the meson-photoproduction processes on a deuteron target provides information about the underlying 
reaction mechanisms on few-body systems. Our project has been motivated by a special interest in the role of 
intermediate particles in reaction (Q. 

In the 1970's, the contribution of intermediate particles and resonances to the differential cross section in the 
backward ird elastic scattering was theoretically discussed in Ref. 1.2). It has been predicted that the contribution 
of intermediate particles, formed in a two-step process, should manifest itself a cusp in the energy dependence of the 
backward differential cross section around the corresponding thresholds. Such an effect, associated with intermediate 
77-meson production, was confirmed by several independent measurements of backward nd elastic scattering |13| . 

The preliminary CLAS photoproduction data 1] give for first time clear evidence for the intermediate rj- meson 
effect. A cusp in the region around the rj photoproduction threshold, E 1 ~ 600 — 800 MeV, is visible in the energy 
dependence of the differential cross section at large CM scattering angles, cos 6* < —0.6. The observed effect becomes 
more pronounced as the scattering angle increases. This behavior was also seen in a previous measurement of reaction 
||TJ 14] in which a small structure was observed in the excitation function the differential cross section at cos 9 — —0.64 
(the maximum scattering angle of this experiment). Reaction (1) was also theoretically considered earlier (see, for 
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I. INTRODUCTION 
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example, Ref. [4j ) in the framework of single- and double-scattering approach without consideration of the intermediate 
r\ meson. A reasonable theoretical description of the previous data |14j was achieved without the necessity of inclusion 
of the u t) effect" at low momentum transfer. 

The aim of the present paper is the theoretical study of reaction at large pion production angles. Our principal 
interest is the contribution of the intermediate r\ meson to the differential cross section and whether it can explain the 
structure in the differential cross section that has been observed in a recent CLAS experiment lj. We shall use the 
standard approach based on single- and double-scattering amplitudes. The main contribution to the total amplitude 
at large angles is expected to come from the double-scattering terms. We shall consider photon energies far from the 
ivd threshold, say 500 < E~ < 1000 MeV, where the influence of the intermediate iV(1535) resonance and the r\ meson 
effect should be important, and it is possible to neglect the excitation of the A(1232) isobar in the intermediate state. 
Thus, we construct the total reaction amplitude from terms, expected to be essential, with approximate values of the 
parameters. In this paper, while we do compare our predictions with the data qualitatively and present all the details 
of our treatment, we do not attempt a detailed description of the CLAS data; we leave this task for forthcoming 
papers from the CLAS Collaboration. 

This paper is organized as follows. In Section [HI we derive the expressions for different terms of the amplitude 
of reaction P in a diagrammatical approach. In Subsection III Al we briefly discuss the main contributions to the 
reaction amplitude and introduce the notation that we use. In Subsections III Bl and III CI we give gauge-invariant 
expressions for the resonance and VME contributions to the elementary amplitudes on the nucleoli. These results are 
used in Subsections III Dl and III El to obtain single- and double-scattering amplitudes of reaction Q, respectively. In 
Section lTTTl we present the numerical results. In Subscction lHI Al we study the influence of the "non-static" corrections 
that are taken into account in the double-scattering amplitude with intermediate r\ production. In Subsection 1111 Bl 
we discuss our numerical results for the differential cross section (its energy behavior at several values of cos 9) of 
reaction (JIJ with backward tt° production. The conclusion is presented in Section IIVI In Appendices lAl iBl [HI 
and [D] we calculate the the integrals and some other expressions, used to derive the single- and double-scattering 
amplitudes of reaction . 

II. FORMALISM 
A. Diagrams and notation 

The diagrams for single- and double-scattering amplitudes Mn ) and M( 2 ) of reaction are shown in Figs. 
and ^3, respectively. The notation for the four- momentum vectors of the initial, intermediate, and final particles are 
given in this figure. The vertices marked by "i " or "j " correspond to the elementary amplitudes of the subreactions 
on the nucleons, and indices "i " , "j " specify the contributions to the elementary amplitudes considered below. In 
Fig. ^3, the notation "/i" stands for the intermediate meson. Hereafter, we shall consider only diagrams with h = tt 
or r\. 

The elementary photoproduction amplitude 7 AT —> hN is usually constructed as a sum of Born, vector-meson- 
exchange (VME), and resonance terms 0, 0, @- The Born amplitudes correspond to a set of tree diagrams with 
NNh coupling and all possible couplings with a photon, summed by a contact 77rAQV-coupling term. It is known 
that the total Born amplitude satisfies gauge invariance (see Ref. Q and references therein) . Using Born amplitudes 
for subprocesses in reaction on the deuteron, one encounters the problem how to get the total gauge-invariant 
amplitude. The problem comes from nucleon off-shell effects, and the way to solve it is discussed, for example, in 
Ref. (see also references therein). However in our analysis, we shall neglect Born terms in the kinematical region 
of reaction under study. 

The resonance and VME terms in thephotoproduction amplitude arc shown graphically in Fig. |3 where VME terms 
are calculated via p and u> exchanges |4j,|9j- The "meson-meson" hN—>TtN amplitudes are written through resonance 
contributions (Fig. |3J) ■ The main contribution from an intermediate r\ meson to the cross section of reaction Q is 
expected from the double-scattering diagram in Fig. ^> (h — 77), with the A(1535) excitation in both blocks of the 
diagram due to the large partial width of the decay iV(1535) — > rjN. Several nucleon resonances |X5i are coupled to 
the ttN system. Also, the couplings of AT(1535), AT(1440) , and A(1520) to the rjN system | HIoTllll | are often used in 
the production amplitudes. However in our qualitative analysis, we shall limit the resonance parts of the elementary 
amplitudes (Figs.|21and EJ) to the contributions of the iV(1535) in the irN and ijN channels and the 2V(1440) in the 
ttN channel. We do not include the contribution of the A(1232) isobar since the considered energies are far above the 
A(1232) region. 

Let us write the total amplitude M 7 <j of reaction JIJ as 
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M ld = M (1) + M (2) , M (1) = M i > M (2) = H M «j . ( 2 ) 

z i,h : j 

where h = ir,rj and index i = 1,2, us, and p (the same for j) stands for iV(1535), iV(1440), uo- and p-exchange for the 
elementary subprocesses, respectively. Note that p-exchange term in the single-scattering amplitude Mm is forbidden 
by isotopic arguments. We shall use standard normalizations for the amplitudes, corresponding to the following 
expression of the differential cross section for a binary reaction: 

-^ = -^W- (3) 



Here, \M\ 2 is the square of the total amplitude, averaged (summed) over the polarizations of the initial (final)particles; 
s is the square of the total CM energy; and q and k are the relative 3-momenta of the initial and final state particles. 
Below, we use the following notation: Lpi 2 (xi,2) are Pauli spinors (isospinors) for the initial and final nucleons in the 
elementary reactions or for the nucleons in the deuteron and tp + ip~ 1 (x + X = 1); o = ( acr )i where a is any 3-vector 
and cr is a 3-vector from the Pauli spin matrices. 

In Subsections III Bl and III CI we give gauge-invariant expressions for the elementary amplitudes on the nucleon, 
notations for coupling constants etc., that will be used for the total amplitude M 7 <2 in Subsections III Dl and III El 

B. Resonance terms in the amplitudes on the nucleon 

In order to obtain the expressions for resonance contributions to the elementary amplitudes, let us use the effective 
Lagrangians for the irNR, r)NR , and jNR interactions: 

LttNR — —ig-irNR NTtR it + h.c. , L V N]i = —ig^NR NTRr) +h.c. , (4) 

L i*R = lu—i rR(k s R + klT 3 )T^NF^ + h.c. (e 2 /A7r « 1/137), (5) 

2(m + m„) 



T= 1, IV=7 5 <7 A „, (odd parity), (6) 
r=7 5 , T I1V = (T IXV (even parity) (7) 

(we use the pseudoscalar couplings in Eq. (0J). Here <j^ v = ^-(7^7^ — 7^7^) ; = d"A^ — d^A" ; N, n, -q, A^, and R 

are the nucleon, 77, ir, photon, and resonance-particle (R) fields; m and m R are the nucleon and resonance masses; k s 
and k v R correspond to isoscalar and isovector jNR couplings. Operator structures in Eqs. © and Q correspond to 

odd- [JV(1535), I{J PC ) = 5(5")] and even- [jV(1440), I(J PC ) = \{\ + )] parity resonances R. 
Hereafter, we use the photon couplings <? 7 j, g S i , and g V i , defined as 

g 7 i = g S i+g V in = — ■ — (ki+kln), (8) 

m+rm 

where i — l[iV(1535)], 2[iV(1440)] specifies resonances and their masses m^. 

Below, we use the non-relativistic deuteron wave function (DWF) to obtain the amplitudes of reaction (JIJ (Sub- 
sections ^TO] and . In this connection, we derive the elementary amplitudes in a non-relativistic approximation 
before they are used to obtain the amplitudes on the deuteron, leaving the leading terms only with respect to the 
relative momenta. The relative accuracy of this approximation near the 77 photoproduction threshold (y/s ~ m+m v ) 
is of the order of q/2m ~ 0.2 — 0.3 (q is the photon relative momentum). Let us introduce some useful notation for 
the resonance amplitude A a n, of the reaction ci-ZV— > R~ > bN as 

Aaib = xt T Xi <p£S<pi , T = T aib , S = Saib, (9) 

where S and T are spin and isospin operators. Let us first consider the spin operators, A^, of the photoproduction 
amplitudes. In our approximation, they have the following structures: 

S 7 ih ~ (?o e - eo if) , S^h ~ k (q e - go eo) , 
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where e = (ere), q = (erg), and k = (crk); k (q) is the CM 3-momentum of the final meson (initial photon); qo and 
e = (eo,e) are CM energy and polarization 4-vector of the photon. Note that gauge invariance, guaranteed by the 
Lagrangian L^nr ©, is obviously valid for the operators S^ih and S 7 2h- Hereafter, we shall fix the photon 4-vector 
by the gauge condition eo = (eq) = 0. Finally, we obtain the following expessions for the resonance amplitudes A a n, in 
terms of the operators T a n, and S a a © : 



(7r 2 T)(7rix) , T ni „ = (tvt) (i = l,2) . 



(10) 



Sjih — iC^ih Fih e , S l2 h — iC 1 2hF2hkqe , Shin — i^hin Fih Fin I C 1 — ^i 7 ?) 



Cyih = gift BWi 2mg , C T 2fe = 32/1 -BPF 2 , = igihgi^BWi 2m, 



S^TT = *Cft27T ^2/l Jp27T ? 7 Cft27T = ~ 02/t 52lT B^l ) BWi — „ j . ^„ , . 

zm s— m\ + ly/s 1 i(s) 



where, in addition to the above-given notations, q is the CM 3-momentum of the initial photon or meson; BWi > r«( s ) , 
and gih are the Breit-Wigner propagator, total width, and coupling constant to the hN channel for i-th resonance, 
respectively; I is the unit 2x2 matrix; and are the form factors of the strong decays Ri — > hN. Here, we use the 
form factor only for p-wave RuNn vertex (Fih = 1 for other vertices) in order to compensate its energy growing in the 
region far away from the irN threshold. We take the function F^n = F in monopole parametrization (see Eqs. iJTTJ) 
which is convenient for analytical calculations of the integrals in Subsection III Dl The widths I\ , coupling constants 
gih , and relative 3-momenta qh (h — 7r, 77) for the decays iV(1535) — > ttN, r/N and iV(1440) — ► ttN are connected by 
the relations 



r x (s) = ri 7r (s)+ri r? (s) , 



_ 2 (E+m) q-x 

J- In — <59ln Z 7= • 
47T JS 



_ 2 (E+m) q v 



(11) 



^2n — 3 5 27r 



J {E - m)Cl * F^), F(qn) = 



A 2 



47T ^/i 



A2 _ a2 , 2 

yv — i\ -rq 0n . 



where E + m 2m, E — m~ q\j(2m) (E is nucleon total energy), and is the relative momentum in the decay 
TV (1440)— >ttN at resonance mass. 

The photon couplings g 7i can be expressed through helicity amplitudes A v x , 2 and A™^ 2 [3 of the decays Ri — > ^7 

and i?i — » tt-7, respectively. We can relate the radiative widths to the amplitudes (for spin-1/2 resonances) as 

well as to the constants <&™ = g S i + g v i, i.e., 



r(iJj-*p7, 717) = 



7 



7T TTti 



\ A P,n |2 _ ,^,» 



where fc 7 is the relative photon 3-momentum in the decay. Then, we obtain 



\A* 



n i 2 f p,n^2 TTlj 

/ 2 I — 



7T7 j /Cr-y 



(12) 



(13) 



Note that only the isovector constants g V { (not g S i) are needed to derive the amplitude for reaction on the 
deuteron. The helicity amplitudes are usually extracted from the photoproduction experiments, and the values A V ^J 2 
for iV(1535), iV(1440) and other nucleon resonances can be found, for example, in Refs. [El llfil|. 

C. Vector-meson exchange (VME) terms 

In order to derive u>- and p-exchange amplitudes of the photoreactions on the nucleon, we use the effective La- 
grangians L VNN and L em of VNN(V — uj,p)- and Vhy (h = n, 7y)-interaction, taken in forms used in Refs. 0,0 



N 



7, + ^<WH<~''' + ^''> 



N . 



(14) 
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(15) 



(£1230 — 1 j e i23 = 1) 5 where p, ^, A, a and i, j are Lorentz and isotopic indices, respectively, and it, ?y, p, and w stand 
for 7r, 77, p, and uj mesons. The coupling constants Gyh-y (h = n, rf) in Eq. Ijl5(l can be expressed through the radiative 
widths IV— >h-y by the following relation: 



r 2 n 3 
\2nml 



G vh 7 ml 



12Trm 2 h 



Tfl 



(16) 



Using Eqs. (fT4*jl and one can write the VME amplitudes of the reactions jN^hN as A 1 vh — xt^Xi U2M 7 vh ui , 
where U1.2 are nucleon Dirac spinors (ui 1 2Ui,2 = 2m), and 



iCr Gyh-f 9vnn a vi \ 



(17) 



where and e v are the 4-momentum and polarization of the photon, k x and p a are 4-momenta of the final meson 
h and the initial nucleon, and r 2 = (q — k) 2 is the 4-momentum transfer to the deuteron squared. T = T 1 vh is the 
isospin operator: 



(ttt) , T 7W7r = (n 3 7r) I, T lpv = T3 . 



(18) 



where 713 = (0,0,1) is the unit vector in isotopic space. Using the notation JHJ for the amplitude A^vh, i-s., 
U2M 7 yhUi — (p^SjVhfi, and the gauge invariance condition eo = (eq) = 0, in non-relativistic approximation, 
we have 



S^ Vh = iC jVh Tr{kqe}I, C jVh = — ^jT 9 ™" , Tr{kqe} = 2i (fe • [gxe]) 



(19) 



(here, the expression for the trace Tr{kqe} is used for convenience in the next Subsection III Djl . The amplitudes 1|19() 
in our approximation do not contain the constants j3 v from Eq. i|14|) ■ The gauge invariance of the amplitude (|17l) is 
obvious and is also satisfied in the expression l|19(l for S^vh due to the gauge-invariant factor Tr{kqe}. 

Note that the pir^ vertex in the Lagrangian 1|15|) corresponds to isoscalar photon coupling, while only isovector 
photon coupling can contribute to the amplitude of reaction Q . Generally, the isotopic piry vertex has the structure 
gi (tvp) + §2 7T3/?3 + 93 («3 [tt x p]), where (73 is an isovector coupling constant. Then, for radiative p-decays, we have 



■*7T U 7 



■ {91+92) 2 and r p ±_ 



g(+gl From the PDG 15], r p o_ w o 7 /r p ±^ 7r±7 ~ 1.7-1.8, i.e., (gi+g2) 2 > 9i+gl 



Since isoscalar pnj coupling (g 2 = 53 = 0) is successfully used in the p-exchange amplitude of the reaction yN-^nN, we 
may suppose that (72.3 <C g\. In addition, let us compare lo- and p-exchange amplitudes M lull and A/ 7p7r of the reaction 
jN^tt°N. Using coupling constants from |4| (gi = G p ^ 7 ), we obtain M lul7r /M. 



G 



~fUJ7\ 
CJ7T7 h)LO 



nn/G 

pnj 9pN N ^ 

10. 

Based on that, we neglect the p-exchange amplitude with intermediate pion production in reaction Q in comparison 
with w-exchange amplitude. 



D. Single-scattering amplitude of the reaction yd — > Tv°d 

Let us write the amplitude A of the process yN — > ir°N in the form 10, where S and T are the spin and isospin 
parts of the transition operator. Then, a single-scattering amplitude for reaction Q reads 

iM (1) = Tr{T} J^--^'TT{i(e^t 2 )iG(p 2 )iSiG( Pl )i(e 1 t 1 )iG(p)} (20) 

(we follow the diagramatical technique of Ref. 0]> and some comments will be given in Subsection III E(l . Here: 
G(pi,2) = (2to£i,2 — p\ 2 + i0)~ 1 and G(p) — (2me — p 2 +z0) _1 are nonrelativistic propagators of the intermediate 
nucleons in Fig.^, with 3-momenta (kinetic energies) p 1 2 and p (£1,2 and s ), respectively; E\ t 2 and T%,2 are 3- vectors 

of polarization and dNN-vertex operators for initial and final deuterons. The vertex T is connected with the DWFs 
by the relations 

(ef) =4V^(g 2 +me d )*, * = a 2V \ xt^xl , (21) 
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where q is the relative 3-momentum of nucleons, ea is the deuteron binding energy, u(q) and w(q) are the s- and 
d-wave parts of DWF normalized as J d 3 q [u 2 (q) + w 2 (q)] — (27r) 3 . 

Integrating Eq. I|20|l over the energy e and using Eqs. (|21|l . we obtain 

M (1) =2Tr{f}y^Tr{i+S* 1 }, *i, 2 = *(ei, 2 , <7i, 2 ) , (22) 

where q 12 = P~\P\,2 (see Fig. |TJi). Then, using Eqs. ©, O, dJl, and (JT^J for the 7iV-^7r°7V amplitudes, we 
obtain 



M (1) = Mi + M 2 + M w , M - x 



Oi=e, <3 2 = fcqe, O w = Tr{fc #e} I , 



Hereafter, go j 9: an d fc are the CM photon energy and the CM 3-momenta of the initial photon and final pion, 
respectively. The values xi j2 , factored out of the integrals in Eq. (|23[) . depend on the effective mass m^N = in the 
subprocess jN irN , and we calculate the value m^jv using the 3-momentum p = i(Pi+P 2 ) of the intermediate 
nuclcon in Fig. [TJi. Expressing DWF, given by Eq. (|21[1 for ^> in r-representation, where 



, fd 3 r _ iar ~, . ~, N u(r)„ w(r) ( 3 (re) „ \ 

(«. = 7 , *( 6 , r) = ^ e - (-^r - c j , (24) 



we obtain the amlitudes Mi (|23|l in the form 



lf,=,,/^ I 6 iA, 'lr{i+d i i 1 }, e il2 = $(e,r), A=i(fc-</i. (25) 



In order to evaluate the amplitudes Mj, let us introduce the integrals 



/ 



" lAr f?(r) = M , e lAr /i, 2 (r)/ 2 (r)^ = n^B^ + S tJ C h2 , (26) 



(4tt) 2 Jiw " J {At:) 2 J ^ ' r 

where 

AW-^ + ^, «r)-M2, o.^(M-l). (27) 
From Eqs. (|23fl - l|27[) . we obtain the amplitudes Mj (i = 1, 2, w) in the form 

Mi =a;i(Ai-2Ci)Tr{e 2 *eei} + a;i-BiTr{Fen}, (28) 

M 2 = x 2 (Ai - 2Ci) Tr{e 2 *fc q e ej + ar 2 -Bi Tr{yfc gen} 

+ x 2 [C 2 (eie 2 ) + {B 2 -2B 1 ) (n Cl )(ne*)] Tr{fcqe} , 

M w = 2x u [(Ai-2Ci+C 2 ) (eiej) + {B 2 -2B 1 ) (n ei )(ne*) ] Tr{fcge} , 

where V— (ne 2 ) ei— (nei) e 2 * (e 2 * = e 2 <x ). Neglecting the d-wave component of DWF, i.e., setting w(r) — in Eq. J57J), 
one obtains i?i j2 = Ci j2 =0 and it simplifies Eqs. (|28|) . However in a single-scattering amplitude, the momentum is 
transferred to one nucleon, and at large angles of the outgoing 7r° the relative momenta q\^ of the nucleons become 
large and the d-wave part of DWF should be important. We use a parametrization of DWF employing the Bonn 
potential (full model) and the corresponding analytical expressions for A\ , Bi 2 , and C\ 2 , used in Eqs. Q28[). are 
given in Appendix lAl 
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E. Double-scattering amplitude of the reaction — > Tv°d 



Let Si(§2) and fi^T^) be spin and isospin operators in the amplitude @ of the subprocess jN^hN (hN^ir°N) 
in the diagram of Fig. Then the double-scattering amplitude Mm) has the form 



zM (2) = Tt{T a T 2 } J — — ^ ^ 

x Tr{i(e 2 * f 2 ) iG(p[) iS x iG{ Pl ) i{e x f x ) iG(p 2 ) iG(p 2 )} iG fc («) . 



(29) 



Here, G(pi. 2 ) and G(pi 2 ) are nonrelativistic propagators of the intermediate nucleons with 3-momenta (kinetic en- 
ergies) p 1 2 and p[ 2 (£1,2 and e[ 2 ); G^(s) is the propagator of the intermediate meson /i with 3-momentum s; 

^2 = t 2T t t 2 , S>2 = (J2S T (J2 19], where index "T" stands for transposition operator. Integrating over the energies e 2 
and and using Eq. (|21|l . we obtain 



(2) 



(2tt) 3 (2tt) 

Inserting T and 5 from Eqs. (JTUJ, l|18|) . and (|19|) . we obtain the contributions Mj^j to M( 2 ) m the form: 

d 3 p 2 d 3 p[ 



(2tt) 3 (2tt) ; 



■■&{0 <w }Gfc(a), 



(30) 



(31) 



Uihi 



2g v 



C 7th C hl7r , Q lhl = tf^etfi , 2 7ri = *Jsge*isfeF(s), 



y ih 2 = —C^ h C h 2^F(k), 6 l7r2 = ^t^iskF(s), 02,2^^isqei> 1 skF 2 { s ) 



2 2 

yuni = G 7W7r G 7r i 7r , y prii — — G 7 p„G„i7r, G W7r i = O p „i = v & 2 Tr{sge} , 

m m 



Vuiir2 



C 1U1 „ F{k) , W7r2 =$+fisfc F(s) Tr{s $ e} , 



where F(s) is a form factor in the iV(1440)7V7r vertex (see Eq. (flip*) . 

In the case of the double-scattering amplitude due to DWF, the main contribution to the integral J d 3 p 2 d :i p' 1 comes 
from the regions p 1 2 ~ hP i and p[ 2 ~ 5-P2 with small relative momenta gi, 2 ~ 0. To simplify the calculations, we 
neglect the d-wave components of DWF in the amplitudes Mihj. The factors in Eqs. (|31|l are calculated at 



Pi=P2 = \ p i, P'i=P2 = \ p 2- 



(32) 



We also factor the value Tr{sge} from oj and p-exchange amplitudes out of integrals (|31|l . fixing the momentum s 
according to prescription (1521 . Then, s = k+ |(Pi— P 2 ) = |(fc+q), and we replace Tr{s ge}-^ |Tr{fc g e} in Eqs. I|31|l. 
The meson propagator Gh{s) — {s 2 ~m\+i0)~ 1 in Eq. (|31|l can be written as 



G h {s) = - 



--(p?+pD 

m 



- iO 



j Qo = <Z0 +Tdl — £d, a h =m h~Qo- 



(33) 



Here, Tdi is the kinetic CM energy of the initial deuteron, s= -^p' 1 —p 2 , and Go is the excess energy in the process 
'yd^pnh. The term (Qo/ TO )(f>i 2 +P 2 ) m Eq. Q33H takes into account the kinetic energies of the intermediate nucleons. 
In this Section, we neglect this term, i.e., we use a "static" approximation. However, in order to study the energy 
dependence of the differential cross section of rj production (Q ~ m„) in the threshold region, we shall restore this 
term in the amplitude with intermediate r\ meson in Subsection IIII Al 
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In the "static" approximation, Gh = — (s 2 + a 2 l —iO) 1 . r-representation is convenient to calculate the double- 
scattering amplitudes. Let us introduce the Fourier transformations for s-dependent parts of the integrands l|3I|) : 

■d 3 v ja _ , ^ sF{s) fd 3 



f err „.„». , „ sF(s) f d A r .„_ T , , , N 
n --e* sr h u (r), 2 W = / — e' ST Vfo 12 r , 34 



-J F%s)= /^ifciW + ^^W] 



S 2 + a 2-i0 4tT ' " 



The functions /in, 12,1, 2(f) are given in Appendix IbI Let us also define the integrals 

J f(r) h n (r) =A U) J f(r) h 12 (r) r = m A 12 , (35) 

/ 2 (r) hi(r)rirj = rrnmj A 221 + fly v4 2 22 , / f 2 {r)r 2 h 2 {r) = A 



1223 j 

where /(r) =u(r)/(rV2), and u(r) is the s-wave part of DWF. For the integrals in Eqs. we use a short-hand 

notation 

(^^-■■ - —I- <»> 

The expressions for j4h , j4i2 , A221 , ^222 , and A223 are given in Appendix [U] Rewriting integrals (|31|l in r- 
representation using Eqs. Q24[l. (|34|) . and (|35|l . we obtain the following expressions for amplitudes M^f 

Mini = yi h i An Tr{e 2 eei} (h = tt,t]), (37) 
Mi n 2 = y\^ 2 A 12 Tr{e 2 *eeimfc} , M 2vl =y 27Tl A 12 TT{e 2 mq eex} , 



M 2w 2 = V2TT2 



A 22 \ Tr{e 2 * mqeeimk} + (A 222 + A 223 ) Tr{e 2 o-qeeicr k} 



Mu, n x= y UJ ^ 1 A 11 (e 1 e 2 )Tr{kqe} , M W3r2 = ^ 2/wtt2 Tr{e 2 * eimfc} Tr{fe qe} , 

M pn i = y prll An (ei e 2 ) Tr{/cq e} . 

By summing the amplitudes from Eqs. (|23|l and l|37|l. we obtain the total amplitude M y d |J2J|. Note that the gauge 
invariance of this amplitude comes from the Lagrangians JHJ and (|15J) and that it is not violated by the nucleon 
off-shell effects in the deuteron. For simplicity, we do not take into account here VME terms in the double scattering 
amplitudes, i.e., the amplitudes M un i , M^ 2 and M pv i will be neglected in the numerical calculations of Section 1TTT1 
The square |M 7< j| 2 , of the total amplitude, averaged (summed) over initial (final) polarizations, is rather cumbersome, 
and we do not write it here. 



III. NUMERICAL RESULTS 
A. Nucleon kinetic energy terms 

Before we consider the differential cross section of reaction JJJ let us discuss the threshold effect of intermediate 77 
production in the double scattering amplitude. The contributions of intermediate r\ meson to the reaction amplitude 
comes from the terms Mi v \ and M m \ l|37|) which contain the integral An- We shall recalculate An in a "non-static" 
case using the r\ propagator l|33(l with nucleon kinetic energy (NKE) terms (Qo/m)(p' 1 2 +p 2 ) taken into account, and 
compare it with the result of the "static" approximation, used in Eqs. I|34l) and l|35|l . 

To simplify the calculation of An for the "non-static" case, let us replace the Bonn DWF by the effective Gaussian 
s-wave function ip(r) = B exp(— 6r 2 ). In order to fix the slope parameter b, let us note that for reaction at the 77 
threshold with backward outgoing ir°, we have a 2 =0 (h = rf) in Eq. (|33|l and Q sw in Eqs. (|35|l . On the other hand, 
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we have An(a 2 =Q = 0) ~<l/r>d from Eqs. 1351) . We fix the value b from the mean value < l/r>d corresponding 
to the s-wave part of the Bonn DWF. The expression of An with the Gaussian DWF is given in Appendix [D] 

Fig. 01 shows Re An and Im^u calculated in the "static" approximation, as a function of the photon laboratory 
energy E 1 for several values of z — cos 9, where 6 is the CM scattering angle of the outgoing ir° meson. Here, one 
can see that the results obtained with Bonn (solid curves) and Gaussian (dashed curves) s DWF are quite close to 
each other. The function Re An peaks at -E 7 = i?fh~630 MeV, where Eth is the threshold for the reaction ^yd^r/pn, 
and ImAn = at E 1 <E t y l . Note that the left (right) derivative d (Re An) j 'dE^ — > +oo [d (Im Aix)/dE 1 — > +oo] at 
E^ — > E t h- These properties come from the "static" approximation used in Eqs. I|34[) . In the r] threshold region, the 

3/2 

function ImAn should depend on E 1 as a 3-particle NNrj phase space, i.e., ~ Q , where Qo is the excess energy 
given in Eqs. l|33|) and Qo ~ Ey — E t h- In fact, when NKE terms {Qo/m){p'i +p\) in the r\ propagator (|3~3l are 

1 /2 

neglected, then ImAn behaves as a 2-particle phase space, i.e., ~ Q Q . 

In Fig. we compare Re An , ImAn , and |An| 2 calculated with the Gaussian DWF in "static" (dashed curves) 
and "non-static" (solid curves) cases. Results are shown for two values z = —0.55 (Fig. |3^,c,e) and z — —0.85 
(Fig. [3),d,f). In the "non-static" case, ImAn ~ Q 3 ^ 2 in the small region close to the r\ threshold and then due to 
DWF begins to decrease where E 1 increases. The energy dependence of |An| 2 clearly demonstrates the ^-threshold 
effect from the loop diagram (Fig. QJ)) in the energy behavior of the differential cross section of reaction when all 
kinematical factors from subreactions on the nucleons are neglected. Fig.|S^,f show that when NKE terms are included 
then | An | 2 turns to be much smoother function instead of exhibiting sharp peaking in the "static" approximation 
case. 

Finally, for the differential cross sections in Subsection IIII Al all double-scattering amplitudes with intermediate 
pion are calculated in a "static" approximation using the Bonn DWF [l^ and the expressions for the integrals 
An l i2,22i,222,223 are given in Appendix IO For the amplitudes with an intermediate r\ meson, NKE terms are taken 
into account and we use An given in Appendix IdI 

B. Differential cross section of the reaction 7c! — > 7T°d 

The amplitude of reaction (JTJ as expressed by Eqs. (|28[1 and l|37|) depends on a number of parameters. In the 
Table [J we list sets of helicity amplitudes A^% for photon couplings to spin-i resonances (see Eqs. (|13|) 'l used in 
our amplitudes. We use the values A^™ from column [20J-1 (1-st variant from Ref. [20J) which is approximately the 

TABLE I: iV(1535) and JV(1440) resonance couplings. Units are (GeV)~ 1/2 x 1(T 3 . 



Ref. 


[20j-l 


[20]-2 


[21] 


m 


[5J 


PDG [15J 


JV(1535) A\ /2 


78 


50 


53 


97 


67 


90 


A" 


-50 


-37 


-98 




-55 


-46 


iV(1440) A p 1/2 


-66 


-64 


-69 




-71 


-65 


A n 
A l/2 


50 


45 


56 




60 


40 



mean values among those given in Table [Q 

For the partial (Tih) and total (Tj) widths l|ll|) of A^(1535) and JV(1440) at nominal masses (y/s — mi. 2), we use 
the values 

Tin =Ti v = 0.5 Ti, r 27r = 0.65r 2 , ri=150MeV, r 2 = 350 MeV, (38) 

and take A= 1 GeV in the form factor F (|llf) of the hadronic iV(1440) decay. 

The coupling constants G UJ1T1 and G prn are obtained through Eq. (|16f) from radiative widths |15|| 

IV^o 7 = 8.7 • 10" 2 T w , r p o_„ 7 = 7.9 • 10~ 4 T p (r u = 8.44 MeV, r p » 150 MeV). (39) 

The strong coupling constants — g^NN and a„ = g„NN are not well determined as it was mentioned in Ref. 0. 
In various analyses, they vary in the ranges [22I l23l| 8 < g u < 20 and 1.8 < g p < 3.2 (we discuss only the vector 
coupling constants, while tensor couplings are neglected in our approximation according to the results of Subsection 
III C|) . Here we list some values of these constants from the above papers: 

2 2 

= 21, g p = 2 [5]; g p = ^=2.66 [4]; g- =0.84, = 20 [18] (full model). (40) 
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For the results, shown below, we take the value g^ from Ref. 01 > *- e -i Sw/4 7r = 20 (g^ ~ 15.85). The constant g p is 
not used in our calculations since p-exchange terms are neglected, as was mentioned in Subsections 111 Al and 111 CI 

In Fig. El we show the calculated differential cross section of reaction Q with backward 7r° photoproduction as 
a function of the photon laboratory energy E 1 at several fixed values of z = cos 9 from z = up to z — —0.85 (the 
experimental CLAS data which are presented in Ref. [lj are at the same values z). The results shown by the solid 
curves obtained with the total amplitude M yd consisting of the terms and (|57|) (without M un i , M U7t2 , and M pr)1 
terms). The other curves are explained in the figure caption. One can see that the contribution from single-scattering 
amplitudes dominates the cross section for z = 0. The relative contribution from the other amplitudes increases as z 
approaches —0.85. Note that the w-exchange amplitude M w dominates in the total contribution from single-scattering 
amplitudes. 

In Fig. |SJ we see a maximum in the energy spectra of the differential cross section at E 1 w 700 MeV, in the angular 
region of z < —0.65. This maximum is getting more pronounced at z— >— 0.85. The CLAS experimental data also 
show the excess of events, but less sharp, of the same order of magnitude (increasing at z— > — 0.85) in a region around 
the same energy for the same angles. For more detailed discussion, we should mention that the effective energies y/s 
in subprocesses on the nucleons in the double-scattering diagram (Fig.HJj) are calculated in the approximation i|32[l . 
These energies y/s are equal to the TyV-threshold mass at E 1 rj 700 MeV. The maxima in Fig. © reflect the rjN- 
threshold effect from the 7V(1535)-propagators in the elementary amplitudes because we use the energy-dependent 
width Ti(s) of V(1535) according to Eqs. {TTJ and Ti v (s) =0 at y/s < m + m v (E.y < 700 MeV). The main effect 
comes from the double-scattering amplitude M\ v i with two V(1535)-propagators while the contribution from Mi^i 
is much smaller due to a small V(1440)V7r-coupling constant (g^ <C gi v )- Thus, Fig. HJ] demonstrates the effects 
of the two-particle (r]N) threshold in the elementary amplitudes on the intermediate nucleons. At the same time 
for the reasons discussed in Subsection IIII Al we do not see a sizable threshold effect from the three-particle (rjNN) 
intermediate state (E t h = 630 MeV). 

Note that the prescription H32[) usually works well due to a rapid momentum dependence of the DWF in comparison 
with ones for the amplitudes of the reactions on the nucleon. However, this approximation does not reproduce 
adequately any sharp peculiarities of elementary amplitudes in the amplitude of a nuclear reaction. Due to "Fermi 
motion" within the deuteron, the r]N threshold in the diagrams in (Fig. 0) are not positioned at a fixed value of 
Ej. These effects tend to be spread over some region of the incoming photon energy. Thus, the sharp maximum at 
E 7 w 700 MeV in Fig. should be smoothed, but we have no proper simple procedure to do this. 

Let us consider the case, when the 7V(1535) width in the elementary amplitudes is a constant nominal value 
Fi(s) = Ti = 150 MeV. In this case, shown in Fig. [7| we have no ?yiV-threshold effects in the elementary amplitudes 
and no sharp maxima at E 1 k, 700 MeV. Here, we see a broad enhancements centered around E 1 — 750 MeV. These 
enhancements appear mainly from the contribution of the amplitude Mi^i with two iV(1535) propagators. Their 
position is shifted to the left of E 1 w 785 MeV (laboratory photon energy on the nucleon target with an effective 
CM energy equal to the JV(1535) mass) due to decreasing factors from the DWF's in reaction QJ. Note that the 
energy-dependent JV(1535) width I\(s) at the rjN threshold (y/s = m + m v ), where r^s) =0, is about two times 
smaller than the nominal value Ti of the V(1535) resonance. Thus, taking iV(1535) with constant nominal width, 
we get smaller values of the corresponding elementary amplitudes in the region close to r/N threshold. Therefore, the 
differential cross sections in Fig. at E~ ~ 700 MeV are essentially enhanced in comparison with Fig. Finally, we 
conclude that the enhancements in the energy region near the r/ threshold, shown in Figs. |H| and [7| are mainly due 
to the V(1535) contributions to the double-scattering diagram with an intermediate rj production. 

The absolute values of the calculated differential cross section at z = (Figs. and |JJ arc in an approximate 
agreement with the CLAS data l]. For larger angles 9, our results are getting lower in absolute values in comparison 
with the data. At z = — 0.85 (Figs.QJ andCf), the calculated differential cross sections are several times smaller than 
the experimental ones. We expect that the contribution from the single-scattering amplitude M u in our treatment is 
too large. On the other hand, the VME terms in the double-scattering amplitudes (not included here), being added, 
may essentially improve our predictions at z — > —0.85, where contribution from the term M u is getting smaller. 
Thus, we hope that decreasing the elementary w-exchange amplitude (taking smaller coupling constant g^NN and 
introducing form factors) and including the double-scattering amplitudes M^-ni and M u »2 (the contribution from 
M prj \ should be much smaller), we may have a good description of the experimental absolute cross sections at z — 
and essentially improve our predictions at z— > —0.85. 

IV. CONCLUSION 

We considered the energy dependence of the differential cross sections of the reaction 7c? — > ir°d in a wide energy 
range around the r\ photoproduction threshold at several backward CM angles of the outgoing pion. Our calculations 
are based on a nonrelativistic diagramatical technique and take into account single- and double-scattering amplitudes. 
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We conclude that the contribution of the double-scattering with intermediate production of an r\ meson can explain 
the broad cusp experimentally observed by the CLAS Collaboration [l| in the energy dependence of the differential 
cross section near the r\ threshold. Indeed, our calculations show that a broad enhancement (with the width of the 
order of 100 McV) appears in the energy behavior of the differential cross section at large scattering angles 9 in the 
^-threshold region. This enhancement becomes more pronounced as 9 increases, and the magnitude of the effect is in 
a qualitative agreement with the CLAS data. 

We discussed the role of the 3-particle TVTV^-threshold effect taking into account "non-static" corrections to the 
77-propagator. Indeed, we found that a sharp energy behavior of the differential cross section at the energy of the 
r) threshold calculated in a "static" approximation was essentially smoothed by taking into account "non-static" 
corrections to the rj propagator. Our calculation results show that the enhancements in the energy dependence of the 
differential cross section are, to a great extent, due to the iV(1535) contributions in the elementary amplitudes of the 
double-scattering diagram with intermediate 77 production. 

Our predictions depend on a number of parameters, some of which are not well established (e.g. the constants 
9pNN and g^NN, form-factor parameters, etc.). Not all of the possible diagrams are considered in our analysis. Our 
calculations do not include VME terms in the double-scattering amplitudes or any other resonance amplitudes besides 
the iV(1535) and iV(1440) contributions. Thus, there is a way for further improvements of our predictions. 
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APPENDIX A: THE INTEGRALS Ai, B 1)2 , AND Ci, 2 , EQS.lt26l) 

We use DWF from Bonn potential 0] given in parametrization of 



m.+q 1 



(A.l) 



,(r) = Yd l e~ m ' r (l + — 



vriiT mfr 2 



where u(q), w(q) (u(r),and w(r)) are s- and d-wave parts of the DWF in q(r)-space (see Eqs. Ij21|l and l|24|l). Using 
a DWF of the type (A.l), one can calculate the integrals in r-space with the help (after integrating over the angle) 
of the formula 



°? n 



exp(ajr) 



r n»+l 



V — V" [S ni - In a t 



n 1 

(th>0), S„ = J]-, ^ = 0, 



(A.2) 



fc=i 



which is valid if this integral converges, i.e., Rea^ < and the integrand is finite at r — > 0. These conditions are 
satisfied our case. 

Using Eqs. I|26() . I|27() . (A.l), and (A.2), we obtain the expression for A\ in the form (in this Appendix for integrals, 
we use a short-hand notation given by Eqs. (A.3)) 



Ji 



Ax = J A 2 (r) = ,h + J 2 + .h (J = J 



d 3 r e , jAr 



u (r) 



E 



Ci cj A 
— -j- arctg 



(4tt) 2 



(A.3) 



(vriij =m t + rrij) , 



u(r)w(r) c id 
V2r 2 



8\/2 7r 



ln(m 2 



3m, 2 - 3A 2 - m? 



• arctg- 



A 



w 2 (r) 
4r 2 



E 

ij 



di dj 3 (m 2 + m) + A 2 ) 2 - 4m 2 m) 
327 



4to 2 to 2 A 

1 J 



arctg- 



rrii 
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For integrals Bi 2 and C1.2 , defined by Eqs. (|26|l and (|27|l . one obtains the relations 

Si+3Ci = / A(r)/ 2 (r) = ^J 2 +3J 3 , Si+Ci= / fi(r)f 2 (r)z 2 - +Jf a , 
B 2 + 3C 2 = f f 2 2 (r)=9J 3 , B 2 + C 2 =Jf(r)z 2 =3K 2 , z = {nr)/r, 



3u(r)w(r) 
2\/2r 2 



Kn = 



3w 2 (r) 



4r 2 



where J 2j 3 are given in Eqs. (A. 3). So, we have 



3 / 1 \ 1/3 

Bi = - f ATi + A 2 - - J 2 - J 3 J , Ci = - f - J 2 +3J 3 -X a -A 2 



5 2 = ^(A 2 -J 3 ), C 2 = I (3J 3 -Jf 2 ), 



Ai = > f^- 



-mj m 



2m 7 A 2 m. 



9n (m 2 -m 2 ) 2 + 3A 4 

- lnfm 2 +A 2 ^ - - — — 

2 \ v' ' 2m 2 A 3 



arete 



A 



k 2 = Y / 3didj 



32 it 



-(m 2 + m 2 )(m 2 -m 2 ) 2 + 3(m 2 + m 2 )A 4 + 2A 6 

4m? m 2 A 3 

1 j 



arete 



A 



(A.4) 



APPENDIX B: THE FUNCTIONS /i n (r), fe 12 (r), fei(r), AND fe(r), EQS. (IMll 



Fourier inversion of Eqs. (|34|l . where form factor -F(s), defined in Eqs. Q1U[) . gives 

-Ar 



/iii(r) = — 
r 

hi(r) =At 
h 2 (r) =-A 4 



, . . .. 2 (l + Ar)e- Ar -(l-ar)e c 
/i 12 (r) = iA 2 ^- 



^3 (A 2 -a 2 ) 



(1 + Ar)e- Ar h (-A,r)-h (a,r) 



(A.5) 



2r 3 (A 2 -a 2 ) (A 2 -a 2 ) 2 

e - Ar (l+Ar)e- Ar -(l-ar)e 
2r 3 (A 2 -a 2 ) + 



fto(a,r) 



3-3ar+a t 



2„2 



r 5 (A 2 -a 2 ) 2 



-a A (a 2 h >0) 
i\a h \ (a 2 h <0) 



APPENDIX C: THE INTEGRALS An , A 12 , A221 , A222 , AND A 223 , EQS (l35ll 

Using Eqs. (A.l), (A.5), and (A. 2), one obtains the expressions for integrals l|35[) . Separately, we express them for 
two cases of parameter a, used in Eqs. (A.5), i.e., a — —a (a 2 >0) and a = ia (a 2 <0), where a=|a/,|. 

1) An ■ For the integral An, we obtain 



2(mij+a) 

Q 



A (Q,a)- L (Q,a) 



(A.6) 



64tt 2 



, . (a+a;)L(a;)+TOijA(x) . . (a+a;)A(x) — rriijL(x) 
o-i(x) = > b\(x) — 



Hereafter in this Appendix, we use the notation 



Ao(x,y) = arctg — , L (x,y) = In [(my +y) 2 + x 2 



rriij+y 



A(x) = arctg , L(x) = - In [m 2 , + (a+x) 2 
ma 2 J 



2) A12 ■ The integral A12 can be written as 



A12 (a) = [ Jl2 ( _ A ) - J 12 («)] , J X 2 (a)=ij f 2 (r) e OT 

(in this Appendix, we use the notation from Eqs. I|35[l and Eqs. (|36[l ). and 



C% Cj 

64tt 2 



'2a 3 -3a 2 m l? +ra 3 , \ , Q , 
^ ^+m 4J rj A (g,a) + |L (Q,a) 



ai(x) = u(x) L(x) + v(x) A(x) , 61 (x) = v(x) L(x) — u(x) A(x) . 



3) A221 and A222 ■ The integrals A221 and A222 can be written as 

-4221 ( a ) = \ {3Z 2 -Zl) , -4 2 22 ( a ) = \ (Zi-Z 2 ) . 



Zi = // 2 (r)r 2 /H(r) 



A* 



A 2 -a 2 



no 



Z11 (— A) — Zn (a) 
A 2 -a 2 



Z w = f(r) 



1+Ar 
2r 



— Ar 



Z u (a) = f(r)r'h (r), 



Z 2 = // 2 (r)(mr) 2 ^(r) 



A* 



A 2 -a 2 



-20 



Z 2 i(-A)-Z 2 i(a) 
A 2 -a 2 



Z20=Jf(r) (mrf 



1+Ar 
2r 3 ' 



-Ar 



^21 (a)= f z (r) (mr)%(r) 



(see /io (r) and /ii (r) in Eqs. (A. 5)). For integrals Z w and Zu one gets 

CiC 



-10 : 



E 



64tt 2 



j 777 - - 

-L (Q,A) + -^A (Q,A) 



Zn(-a)=Y, 



Ci Cj 

64tt 2 



Zn (ia)=^^ [oi(Q) + 01 (-Q) + ih(Q) + *6i(-Q)] 



64tt 2 
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ai(x) — u(x) L(x) + v(x) A(x) , b\(x) = v(x) L(x) — u(x) A(x) , 



u(x) = i (x 2 -a 2 -3m ij ) , v(x) = ^ (3x 2 - a 2 - m^) . 



For integrals Z 2 o and Z 2 i, we have 
Z 20 = 5Z 



64tt 2 



""•< ± io( g,A) + ^- 2A)(m ^ +A)2 



2niimj 1 
3g 2 ~ 6 



3Q 3 



4)(Q,A) 



(A12) 



64tt 2 



m i m j f 7a 2 i Smjjd^ i f 3Q 2 a 2 — 3nii 



Q 2 V 15 



10 



10 



£o(0) 



'3m, 3Q + m i3 - (3r4 - 10m?- o 2 + 15a 4 ) + 8a 5 \ 



Z 21 (m)=£ 



64tt 2 



30 Q 3 

ra,ra, / 7 a 2 .3mj,a 



o 2 



15 



10 



+ ai(Q) + ai(-Q) + ibi(Q) + ifti(-Q) 



ai(x) = i(a;) + u(x) A(x) , 6i(a;) = w(a;) L(x) — u(x) A{x) , 



u(x) 



3x 2 a 2 + 3 to 2 , 2a 5 
'- + - 



10 6 ' 15x 3 

4) A 2 23- The integral A 223 can be written as 



v(x) = ^x+^f [3m 4 + 10m 2 - 2 + 15a 4 ] 
4 oO x 



A 22S (a)= Jf 2 (r)r 2 h 2 (r)= J ^ 



J30 + 



J 3 i(-A)-J 3 i(a) 



A 2 -a 2 



(A13) 



30 



/ 2 (r) 



-Ar 



2r 



Jsi(a)= //V) 



1 — ar 



(see /12 (r) in Eqs. (A. 5)). The expressions for J 30 and J31 have the forms 

2(m y +A) 



30 



1 Ci Cj 

2 ^ 64tt 2 

» j 



Q 



4>(Q,A)-L (Q,A) 



J 31 (-a)=£ 



64tt 2 



+ ^ Uo(Q,a) + f m«Q+ ^ * U (Q,a) 



J 31 ( io) = ^ [oi(Q) + ai(-Q) + 161 (Q) + M-Q)} , 



(A14) 



ai (a;) = u(x) L(x) + v(x) A{x) , b\ (x) = v(x) L(x) — u(x) A{x) , 
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APPENDIX D: THE INTEGRAL An WITH NKE TERMS 



The integral An can be written in p-space as 

d 3 p 2 d 3 p[ ip(qi)<p(q 2 ) 



An 



(2tt) 3 (2tt) 3 2D 



, D = (p , 1 +p 2 y+ 1 (p 2 1 +p 2 2 ) + al-iO, 7 



Qo 



(A15) 



(see a 2 and Qo m Eqs. (J22J)), where p' x , p 2 are CM 3-momenta of intermediate nucleons in Fig. [T]p, q x —p 2 + ^q, 

i 



q 2 = —p' 1 — ik and ^(91,2) are s DWF's. We use the Gaussian DWF given in p- and r-space by the relations: 



tp(p) = A exp(-ap ) , ip(r)=B exp(-6r ) , b- 



4a 



(A16) 



Air (rrn+mj) ' 



where N$ is normalization of s-wave part of the Bonn DWF, and c, and mi are the parameters from Eqs. (A.l). For 
denominator D in Eqs. (A. 15), we can use representation 



' ' - iDt dt. 



D =lle 



Then, calculating the Gaussian integrals Jd 3 p' 1 d 3 p 2 , we finally obtain 

00 

f iA 2 
A u = / dt J(t) , J(t) = _°_„ /9 exp 



2 (4tt) 3 D 3 /2 



iatCo o 

ia h t 

ADn h 



(AA7) 



D Q = a 2 - 7 (2+ 7 ) t 2 + 2ia (l+ 7 ) t, C = 4a Q 2 + 7 [a (g 2 + fc 2 ) + it (2+7) fc 2 ] , 

where Q — ^\q+k\. When 0^ = ("static" approximation) and Q = 0, one can express An through <l/r>d- Writing 
An in r-space, we obtain 



Nc /I 



(AA8) 



We use the Bonn DWF in order to obtain the value < 1/r >d and get the slope parameter b of the effective Gaussian 
DWF (A. 16) through the relations: 
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a) Mi b) M ihj 

FIG. 1: Feynman diagrams for the jd^iv d reaction, (a) single- and (b) double-scattering. 





FIG. 2: Diagrams of meson photoproduction on the nucleon (resonance and VME contributions). 




FIG. 3: Diagrams of meson-nucleon binary reactions (resonance contributions). 
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FIG. 4: (a), (c), and (e) [(b), (d), and (f)] show real [image] parts of the loop integral An of the double-scattering diagram 
(Fig. ^j) in "static" approximation with intermediate rj meson in reaction The results are given for the values (a) and 
(b) z — cos 6 = 0; (c) and (d) z — —0.55; and (e) and (f) z = —0.85. Solid (dashed) curves correspond to the results with the 
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FIG. 7: The same as in Fig. [§]but the results are obtained with constant total width of iV(1535). 



